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Abstract 

We present a time-dcpcndcnt Dalitz plot measurement of CP violation parameters in — > 
KgTr~^Tr~ decays. These results are obtained from a large data sample that contains 657 x 10® BB 
pairs collected at the T{4S) resonance with the Belle detector at the KEKB asymmetric-energy 
e+e~ collider. For the CP violation parameters, we obtain two consistent solutions that describe 
the data well. The first of these solutions may be preferred by external information from other 
measurements. There is no evidence for direct CP violation in B° ^ /9°(770)A'g, B° /o(980)i^g 
and B^ i^*+(892)7r^, while measurements of mixing-induced CP violation in B^ p^{770)K^ 
and B^ — > fo{980)Kg decays are consistent with that of 6 — ccs decays. We also measured the 
phase difference between S° — if*+(892)7r~ and 5° — K*~(892)7r+, which may be used to extract 

03- 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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INTRODUCTION 



CP violation in the Standard Model (SM) is due to a complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [H, 0]- At present, mixing-induced CP 
violation has been clearly observed by the BaBar and Belle [4] Collaborations in the h 
CCS induced decay, ^ J/if) Kg, while many other modes provide additional information 
on CP violating parameters. Of recent interest is CP violation in 6 ^ qqs transitions, which 
proceeds by loop diagrams that may be affected by new particles in various extensions of 
the SM. Furthermore, the CP asymmetries in 6 — > qqs transitions are expected in the SM 
to be slightly higher than those observed in 6 — ccs transitions jl, 0, 0, B 1, 10, 11, 12, 13 . 
However, current experimental measurements pW^] tend to be lower than those for b — >■ ccs 
transitions motivating more precise experimental determinations. 

This analysis is concerned with time-dependent CP violation in — > Kgir^n^ decays 



including the quasi-two-body modes P° — * p^{770)Kg and B^ 



s 

/o(980)K^. Quasi-two- 
body time- dependent analyses have been performed on P" p^{770)Kg by the BaBar 
Collaboration [l5| and P° /o(980)irg by both the BaBar and Belle Collaborations [l6|,|l7 . 
However, a two-body approach to these modes is not ideal due to interference with other 
resonances as well as non-resonant decays into the same final state. These effects can be 
taken into account with a time-dependent Dalitz plot fit. In addition, the interference terms 
in some cases may be sensitive to the cosine of the effective weak phase difference {(f)f^) 
in mixing, potentially resolving the two-fold ambiguity in (pf^ characteristic of quasi-two- 
body analyses. Current direct measurements of cos 2(/)i from b — > ccs prefer the solution of 
21.5° over 



18, nj, 20, 21 



22l . |23| | and similarly, it is important to attempt 
qqs modes as well. 

TT- decays is P° ^ K*+{892)n-. Since 



Kin 



/l — uvci yji — 68.5° 

to resolve this quadratic ambiguity in b 

Another mode of interest included in P° 
direct CP violation has been observed in B^ K^n^ (24, l25|, CP violation could also 
be present in B^ i^*+(892)7r^ decays, though existing measurements are consistent with 
zero ' 
in P° 



26| 



This measurement may also be related to the AcpiKir) puzzle; the asymmetry 
K~^7r~ is found to be significantly different from that in B^ — >■ K^ti'^. This 
unexpected result may indicate the presence of new physics (NP) or poor understanding 
of strong interaction effects in B decays. A model-independent test for NP is possible via 
an isospin sum rule with high statistics [13]. Furthermore, it has been proposed that the 
phase difference between P° K*~^{892)7t~ and P° — >■ ii'*~(892)7r"'" can be used to extract 
The BaBar collaboration has also performed a time -dep endent Dalitz 
TT^ decays and released a preliminary result [32 



28, 29, m 31 



plot analysis on P° 



0^+ 



DECAY AMPLITUDE 

The Dalitz plot variables are defined as the invariant squared masses, 

s± = ip±+PoY, (1) 

where p± (po) sue the 4-momenta of the vr^ (Kg) respectively. The final combination, sq = 
{p+ can be obtained from 4-vector conservation, 

■So = ""^lo + 2m^+ + m^o - s+ - s_. (2) 
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The differential decay width with respect to the Dahtz plot variables is 



1 I/4P 

dr(B"^A-V.-)=p-pXl^dM._. (3) 

where A is the Lorentz-invariant amplitude of the decay. In the isobar approximation, the 
total amplitude of B^[B^) K^tt~^'k~ is written as the sum of intermediate decay channel 
amplitudes with the same final state, 

A{s+,s_) = Y.a[Fi{s+,s_), A{s_,s+) = J2a[Fi{s^,s+), (4) 

i i 

where a'^ = \a[\e^^'' are complex coefficients describing the relative magnitudes and phases 
between the decay channels, which also carry the weak phase dependence. The amplitudes, 
contain only strong dynamics and thus, Fj(s_|_,s_) = They can be 

expanded in terms of invariant mass and angular distribution probabilities, 

Fl^is^,s.)^Xtip*) X X Z^ip,q) x (5) 

where p* is the momentum of the bachelor particle in the B^ rest frame, p and q are the 
momenta of the bachelor particle and one of the resonance daughters in the resonance's 
rest frame respectively, L is the orbital angular momentum between the resonance and the 
bachelor particle while are the Blatt-Weisskopf barrier factors fss']. As these factors 
are unknown in general and the Blatt-Weisskopf forms do not improve this model, they are 
taken to be unity with the uncertainty in this choice being treated in the systematic errors. 
The angular distribution, Z^ijI^q) depends on L, 

Zl{p,q)=-^P- g, 

Z'Mq)=lW-q)'-{\mn (6) 

The mass shapes are denoted as -Rj(s+, s_), which differ depending on the decay channel. We 
utilize the Relativistic Breit-Wigner (RBW) js^], Gounaris-Sakurai (GS) [35|] and Flatte [Si^ 
line shapes. Table |T] summarizes the components considered in the signal model, which was 
motivated by the previous Belle Dalitz plot analysis of B^ —>■ Kgn^-n-^ [37]. The fv (1300) 
resonance of unknown spin that appears in the table was first described in Ref. |37| and is 
assumed to be a scalar. 



TIME-DEPENDENCE 

The decay of the T(45') produces a BB pair of which one (B^^^) may be reconstructed 
as B^ —>■ KgTT^TT^ while the other (-Bxag) ^lay reveal its flavor. The proper time interval 
between the -B^ec -^Tag mesons that decay at times, tRec and txag, respectively, is defined 
as At = t-Rec — ^Tag- For coherent BB production in a T(4S') decay, the time- dependent 
decay rate when -Bxag possesses flavor, q {B^: q = +1, B'^: g = — 1), is given by 

+ - q{\A\^ - \A\^) cos Am^At 

(7) 



|A(At,g)p 



4rR0 



+2q^{AA*) sin Am^At 



5 



TABLE I: Summary of the resonances considered in the signal model. All fixed parameters are 



taken from Ref. 
from Refs. 



3411 w ith the exception of those for the /o(980) and /x(1300), which are taken 



38l | and [33], respectively. 



Resonance 



K*±(892) 
7^0^(1430) 
p0(770) 
/o(980) 

/2(1270) 

/x(1300) 

(Kg7r+)NR7r- 
(Kg7r-)NR7r+ 
(vr+7r-)NRi^g 



Fixed Parameters (MeV/c^ 



Form Factor, i?j(s+,s_ 



m = 891.7 ±0.3 

r = 50.8 ±0.9 

m= 1414 ±6 

r = 290 ± 21 
m = 775.5 ± 0.3 
r = 146.2 ±0.7 

m = 965 ± 10 
= 0.165 ±0.018 GeV'^/c'^ 
gK = (4.21 ± 0.09)5. GeVVc^ 

m = 1275 ± 1 

r = 185 ± 3 
m = 1449 ± 13 

r = 126 ± 25 



RBW 
RBW 

GS 
Flatte 

RBW 
RBW 



where r^o is the lifetime and Am^ is the mass difference between the two mass eigenstates 
of the neutral B meson. This assumes no CP violation in mixing, \q/p\ = 1, and that the 
total decay rate difference between the two mass eigenstates is negligible. The amplitudes, 
A, were defined previously; we choose a convention where the B^-B^ mixing phase of g/p is 
absorbed into the B^ decay amplitude, a'^. 

These complex coefficients can be redefined in a way that depends on the decay amplitude, 

a[ = ai{l + Ci)e'^^^+'^^'> (8) 

for A and, 

a: = a,(l - Q)e^(''»-'^») (9) 
for A, and thus a resonance, i, has a direct CP violation asymmetry given by 



For a CP eigenstate, the CKM angle, 4>f^{i), is reduced to a fit parameter, 

- = (11) 



and its effective mixing-induced CP violation asymmetry is calculated as 

|<|2±|a^|2 l + c: 



- n^^cpi^) = iLm = sin20f (z), (12) 
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where r]i is the CP eigenvalue of the final state. Note that Acpii) and S(^p{i) are restricted 
by these definitions to lie in the physical region. For fiavor-specific states the phase difference 
is calculated as 

A(j){i) = arg^af) = 2d,. (13) 

DATA SET AND BELLE DETECTOR 

This time-dependent Dalitz plot measurement of CP violating parameters in 5° 
Kg7r~^7r~ is based on a data sample that contains 657 x 10^ BB pairs collected with the 
Belle detector at the KEKB asymmetric-energy e'^e^ (3.5 on 8 GeV) collider [39|. Operat- 
ing with a peak luminosity that exceeds 1.7 x 10^"^ cm~^s~^, the collider produces the T(4S') 
resonance {^/s = 10.58 GeV) with a Lorentz boost of = 0.425, opposite to the positron 
beam direction, z, which usually decays into a BB pair. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of Csl (Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere jiO]. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector (SVDl) was used for 
the first sample of 152 xlO^ BB pairs, while a 1.5 cm beampipe, a 4-layer silicon detector 
(SVD2) and a small-cell inner drift chamber was used to record the remaining 505 xlO^ BB 
pairs [4l[. We use a GEANT-based Monte Carlo (MC) simulation to model the response of 



the detector and determine its acceptance [42 



EVENT SELECTION 

We reconstruct B candidates from a candidate and a pair of oppositely-charged tracks. 
Charged tracks satisfy loose criteria on their impact parameters relative to the interaction 
point (IP), dr < 0.4 cm and \dz\ < 5.0 cm, where r is the radial coordinate of the Belle 
detector. With information obtained from the CDC, ACC and TOF, particle identification 
(PID) is determined with the likelihood ratio, A/(A + Cj). Here, £j (£j) is the likelihood 
that the particle is of type i (j). To suppress background from particle misidentification, 
vetoes are applied on particles consistent with kaon, electron or proton hypotheses. The 
transverse momenta of charged tracks are required to be greater than 100 MeV and the 



additional SVD requirements of two z hits and one r — cj) hit [43| are imposed so that a 
good quality vertex of the reconstructed B candidate can be determined. We only consider 
Kg —>■ 7T~^7i~ candidates with vertices that are displaced from the IP and that lie in the mass 
window, |m(7r+7r") — m{Kg)\ < 15 MeV/c^. 

Reconstructed B candidates are described with two kinematic variables: the beam- 



constrained mass, Mbc = \/(-E'bc^m)^ ~ (p^^^)^ and the energy difference, AE = E'^^^ — 
-^beam whcrc -Ebeam the beam energy and E^^^ {p'b^^) is the energy (momentum) of the 
B meson all evaluated in the center-of-mass system (CMS). The signal region is defined as, 
5.27 GeV/c2 < Mbc < 5.29 GeV/c^ and -0.04 GeV < AE < 0.04 GeV. The dominant 
background in the reconstruction of -B^ec is from continuum {e'^e~' —* qq) events. Since 
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their topology tends to be jet-like in contrast to the spherical BB decay, continuum can 



be suppressed with a Fisher discriminant based on modified Fox- Wolfram moments [44 
This discriminant is combined with the polar angle of the B candidate in the CMS, cos6b, 
which follows a 1 — cos^ 6b distribution for BB events while being fiat for continuum. A 
requirement that rejects 91% of the continuum background while retaining 73% of the signal 
events is applied. The second largest background comes from charm decays of the B meson 
and peaks in the signal region. We apply the charm vetoes summarized in Table [ITl The 
fraction of events having more than one reconstructed candidate in regions dominanted by 
particular quasi-two-body modes is 28% for B° p°(770)i^g, 25% for B° /o(980)ir^ and 
21% for B^ —>■ K*^{892)7T^. Selecting the B candidate having an Mbc value closest to the 
nominal B meson mass, the fraction of misreconstructed events is 3% for B^ — ^ p^[770)Kg, 
2% for B° /o(980)i^g and 2% for B° K*+{892)7i-. This criteria introduces negligible 
bias into the AE distribution. From this best candidate, the Dalitz plot coordinates in the 
signal region, s±, are calculated after a mass-constrained fit of the B meson to ensure all 
the events are within the Dalitz plot boundaries. 

TABLE II: Summary of charm vetoes. The subscript in the region vetoed indicates that an alternate 
mass hypothesis has been applied to the pion candidates used to calculate the invariant mass term. 

Region vetoed Mode vetoed 

\m{K'^TT^) - m{D+)\ < 100 MeV/c^ B° L'+[A'g7r+] vr" 

\m{K^Tr^)K - m{D+)\ < 15 MeV/c^ B° D+[K'^K+] vr" 
|m(7r+^-)^ - m{J/ip)\ < 70 MeV/c^ J/^P[n+^^-] 

|m(7r+7r-) - m(xco)| < 70 MeV/c^ 5° ^ Xco[vr+^"] 

|m(7r+7r-)^ - m{^{2S))\ < 50 MeV/c^ 5° ^ ij{2S)[n+ fi'] 

Since the B^^^ and B^ mesons are approximately at rest in the T(45') CMS, the differ- 
ence in decay time between the two B mesons. At, can be determined from the displacement 
in z between the final state decay vertices. 

At ~ (^R'^c - ^Tag) ^ 

The vertex of reconstructed B candidates is determined from the charged daughters using 
the known IP. The IP profile is smeared in the plane perpendicular to z to account for the 
finite flight length of the B meson in this plane. To obtain the At distribution, we recon- 
struct the tag-side vertex from the tracks not used to reconstruct B^^^ j43| and employ the 
flavor tagging routine described in Ref . [45| . The tagging information is represented by two 
parameters, the -Bxag flavor, q and r. The parameter r is an event-by-event, MC determined 
flavor-tagging dilution factor that ranges from r = for no flavor discrimination to r = 1 
for unambiguous flavor assignment. The total effective tagging efficiency is determined to 
be 0.29 ±0.01. 



SIGNAL YIELD 

The signal yield is extracted with an extended one-dimensional unbinned maximum like- 
hhood fit to AE with the requirement, 5.27 GeV/c^ < Mbc < 5.29 GeV/c^. The signal 
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FIG. 1: Fit result for — > K^vr+vr signal yield. The solid curve shows the fit result, the dotted 
curve represents BB and the dashed curve is the sum of the qq and BB components. 



probability density function (PDF), "Psig, is modeled with a sum of two Gaussians, where 
the fraction and shape parameters of the tail Gaussian are fixed from MC relative to the 
main Gaussian. The qq component, Vgq, is modeled with a first-order polynomial whose 
slope is allowed to float in the fit while the BB component, V^b is described by a histogram 
determined from MC. The total likelihood for events becomes 

= 77^ n Nsi,Vsi,{AE,) + N^^,[f,,V,,{AE,) + (1 - /,,)P^5(AE,)], (15) 

■ i=i 

where the qq fraction, fgg, is a free parameter. A fit to data yields Nsig = 1944 ± 98 
events with the fit result shown in Fig. [TJ In the signal region, there are a total of 4547 
5° — i> Kgir'^n' candidates that are used to perform the time-dependent Dalitz plot analysis. 
The average signal purity in this region is 0.41 ± 0.02 and the contribution of continuum to 
the total background is found to be 0.89 ± 0.02. 



SQUARE DALITZ PLOT 

In 5° — i> KgTT^Ti^ decays, both signal and continuum events tend to concentrate near 
the Dalitz plot kinematic boundaries making large variations in these small areas difficult 
to describe if non-parametric shapes such as histograms are used to describe the continuum 
background. In order to alleviate this problem, a transformation that expands these regions, 

ds+ds^ — > \det J \dm'd9' , (16) 

is introduced to define the Square Dalitz Plot [46]. Its coordinates are 

m = — arccos 2 —. 1 , = —On, (17) 
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where mo is the vr+vr invariant mass, m™™ = 2m.n 



itlbo — m^o are the kinematic 



hmits and 6*0 is the angle between the vr and in the vt+tt frame. The determinant of 
the Jacobian of this transformation is 



|detJ| = 4|p_||po|"^o 



sin vrm' sin nO' , 



(18) 



with \p_ 



El — and |po| 



El — m'^j^o evaluated in the tt+tt frame. 



SIGNAL PDF 



The signal PDF is given by 
^True(m', 0', At, q) = e(m', ^')|det J| 



^-\At\/TgO 



{I - q/\w){\A\' + \A\') 

- 2w){\A\^ - lAp) cos Am^At 
+2g(l - 2w)^{AA*) sin Am^At O R 



BO BO, 



(19) 



which accounts for CP dilution from the probability of incorrect flavor tagging, w and 
the wrong tag difference between and B^, Aw, both of which are determined from 
flavor specific control samples 45|]. This PDF is convolved with the At resolution function, 
R^obo, described in Ref. ji^]. On the other hand, detector resolution is ignored in the 
Dalitz plot because the widths of the dominant resonances are large compared to the mass 
resolution. For example, the detector resolution in the region of the narrowest resonance, 
the /o(980), is typically 4 MeV, while most estimates of the /o(980) width tend to be at 



least an order of magnitude larger [3J]. The relative signal detection efficiency variation 
across the Dalitz plane due to detector acceptance, e(m', 9'), is also considered and is taken 
from MC. Misreconstructed events usually occur around the corners of the Dalitz plot with 
a slow pion from the tag side. The fraction of misreconstruction can be as high as 10% in 
these regions. As misreconstruction is not accounted for in the signal model, its effect will 
be included in the systematic errors. 



CONTINUUM PDF 

The PDF for the time-dependent Dalitz plot for the qq continuum is determined from a 
narrow sideband, 5.250 GeV/c^ < Mbc < 5.265 GeV/c^ and -0.1 GeV < AE < 0.1 GeV, 
which ensures that the resonances present in continuum are not shifted significantly as the 
Dalitz variables are scaled to originate from B decays. There is a flavor tag asymmetry in 6' 
due to the jet-hke topology of continuum; as a high momentum 7i~^{7i~) in B^^^ is correlated 
with a high momentum 7r^(7r+) on the tag side. To account for this, we fit the Dalitz plot 
asymmetry with an empirical two-dimensional PDF, 

Agg{m', 9'; B, C) = [m\B + Cm')]{9' - 0.5). (20) 

The Dalitz plot component of the qq PDF will be modeled with a two-dimensional histogram, 
Hqqim' , 9'), including a term accounting for the flavor tag asymmetry with parameters fixed 
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from the two-dimensional fit to sideband. To reduce the effects of statistical fluctuations 
with limited data, the symmetry condition, Hqq{rri\ 9') = Hqq{m', 1 — 6') is imposed. 
The At distribution is modeled with a lifetime and prompt component, 

-|At|/r^^- 

Pqq{At) ^ (1 - /,)— + n S{At - fis), (21) 

convolved with a double Gaussian, 

Rqq{At) = (1 - ft,n)G{At 

A^mean) 'S'main'S'tail)) (^2) 

which uses the event-dependent At error constructed from the vertex resolution, a = 
(\/'^Rec + '^Tag) / /^T^ ^ scale factor. The total qq PDF is given by 

Vqq{m', 0', At, q) = l + e')PqqiAt). (23) 

As a check of this PDF, we perform fits in various Dalitz plot regions to extract the At 
shape, however no correlation between At and the Dalitz plot coordinates was found. 

CHARGED B PDF 

Like continuum, B^B~ events also exhibit a flavor tag asymmetry. Using a generic MC 
sample that contains charm and rare B decays, a similar procedure is performed to model 
the asymmetry in the Dalitz plane. The Dalitz component of the B^B" PDF is modeled 
with a two-dimensional histogram symmetrized as Hs+B-ijn' i = i7s+B-(m', 1 — 6''). 

Since charged B events are misreconstructed by borrowing a particle from the tag side, 
the average At lifetime tends to be smaller and should be taken into account. The At 
distribution is modeled with an exponential PDF with effective lifetimes determined from 
MC, convolved with the At resolution function for B^, Rb+b-, 

g-|At|/Teff 

Pb+b- (At) = — ® Rb+b- ■ (24) 

^Teff 

The B^B~ time-dependent Dalitz plot PDF becomes 

VB+B-{m',e',At,q) = i±^::^^^l^^i/s+B-(m',^')^'B+B-(At). (25) 

NEUTRAL B PDF 

As in the charged B case, this PDF is also fixed from MC. However, no correlation between 
flavor and the Dalitz coordinate is found. A two-dimensional histogram is chosen to represent 
the Dalitz distribution that again incorporates the symmetry condition, i7goso(m', ^') = 

i/Bo^o(m',l-^')- 

Lifetime fits to the At distribution are also consistent with the nominal lifetime value, 
therefore the current world-average B^ lifetime will be used as the default value in all fits. 
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g-IAtl/r^O 

4:TbO 



As neutral BB pairs mix, a time-dependent CP violating PDF is used to model At and q, 

1 — qAw + g(l — 2w) (^Abo bo cos Am^At + Sbo bo sin AmdAt^ (g) Rbobo, 

(26) 

with ^BOfiO = ^BOBO = 0. The total B^B^ PDF is, 

PsofioK, e', At, q) = HBOBo{m, e')PBOBo{At, q). (27) 

Significant B^B^ backgrounds with measured CP parameters are treated separately using 
Eq. EH There are ~ 20 5° ^ r/'i^^ and ~ 4 ^ ai(1260)+7r- events expected in the 
data sample. In these decay channels no flavor tag asymmetry in MC is found and lifetime 
measurements are consistent with the nominal value. The CP parameters for B*^ — > rj'K^ 



and B^ — > ai(1260)"*"7r are taken from Ref. 34 



OUTLIER PDF 

To account for events with At values not yet described by either signal or background 
PDFs, an outlier PDF is introduced, 

Pout(m', e\ At, q) = ^H{m', e')C{At; 0, dout), (28) 



where the width of the Gaussian, aout; is determined in Ref. [43| and H{m' , 6') is the Dalitz 
plot histogram of data. 



TIME-DEPENDENT DALITZ PLOT PDF 

The full time-dependent Dalitz plot PDF is given by 

V{m,e',At,q; AE,r) = 

(1 - /out){/sig(AE,r)Psig(m',^', At,g) + 

[1 - /sig(AE, r)][f^,{AE)V,,{m', 6', At, q) + 

{1 - U,{AE)]{fB+B-VB+B- {m', 6', At, q) + 
fBOBoVBOBoijn' , e'. At, q) + 
fv'KoV^'Ko{m',0',At,q) + 

(1 - Ib+B- - Ibobo - frj'Ro) 
xVa,(1260)+n-{m',0',At,q))]} + 

foutVout{m',9',At,q), (29) 

where /gig(A£',r) is the event-dependent signal probability, 

f ^ A P ^) ^ Psig(r)Psig(AE) . 

^'''^ ' ^ psig(r)Psig(Ai?) + [l-psig(r)]K-n,(AE) + (l-p,,)P^5(Ai?)]- ^ ' 
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The r-bin dependent purity, psig{r), is calculated in the signal region and is included to 
increase signal sensitivity and Pgg = Nqq/{Nqq + A^^^) is the qq contribution to the total 
background in the signal region. The event- dependent qq probability is calculated as 

^''^^^^ ^ pMAl)lt-pl)VsB{^Ey ^^^^ 

and /b+b-, Ib^b'^i fr)'K° are constants defined in the signal region as 

N- 

f. = _1 ^ (32) 

Nb+B- + N^OBO + N^,j^O + A^ai(1260) + 7r- ' 

where A^, are the expected number of events in background category i. 

As there is only sensitivity to the relative amplitudes and phases between decay modes, 
we fix ax*+(892) = 1 and 6x*+(892) = 0. In addition, the /2(1270), /x(1300) and non-resonant 
components share common CP parameters giving a total of 27 free parameters. 



FIT RESULT 

We perform a time-dependent Dalitz plot fit to data and find four solutions with consistent 
CP parameters given in Table IIIII These were obtained by performing numerous fits with 
randomly generated initial values for the free parameters. For each resonance, i, a relative 
fraction can be calculated as 

I{\A\^ + \A\^)ds+ds^ ' ^ ^ 

where the sum of fractions over all decay channels may not be 100% due to interference. 
Table HVl summarizes the relative fractions for all solutions. 

Solutions 1 and 2 were found in the previous Belle analysis ^3] and correspond to the two 
solutions found in the higher statistics Belle Dalitz plot analysis of the isospin partner mode, 
K^-K^-K^ [47]. The signal models used in this analysis and that of the isospin partner 
are almost identical. These two solutions are due to an interplay between the two broad 
S-wave amplitudes, /ro^(1430)7r~ and the non-resonant component, and are characterized 
by their different relative fractions. However, they give almost indistinguishable Dalitz plot 
distributions; studies with high statistics pseudo-experiments confirm the existence of these 
two solutions. Solution 3 has a high /o(980)i^'5 fraction compared with Solution 2 that may 
indicate a similar type of interference between /o(980) and the non-resonant component. 
Unlike the previous case. Solution 3 cannot be found in a high statistics toy MC sample 
generated with Solution 2, neither can Solution 2 be reproduced from Solution 3 indicating 
that one of these solutions cannot survive with higher experimental statistics. Considering 
that a solution with a high /o(980)i^5 fraction was not found in the previous analysis or 
the higher statistics K^7r^n~ analysis [47], we therefore conclude that Solution 3 is 

due to statistical fiuctuations in the small B^ — > iiT^vr+Tr" sample and is unlikely to be a 
physical solution. Solution 4 appears to be statistically consistent with Solution 2 for all 
parameters, while Solution 4 cannot be reproduced from a high statistics toy MC sample 
generated with Solution 2 and vice versa. This indicates that these two solutions are likely 
due to a statistical fiuctuation of a single solution and we select Solution 2 as the nominal 
solution of the two because of its slightly better likelihood. 
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TABLE III: Time-dependent Dalitz plot fit result. The phases, 6^ and di, are given in radians. 



neter 


ion 1 


ion 2 


ion 3 


ion 4 


Parai 












'o 

CO 


en 




Cx*+(892) 


0.106 ± 0.058 


0.099 ± 0.057 


0.081 ± 0.058 


0.116 ±0.058 


dK*+ (892) 


-0.006 ± 0.200 


0.127 ±0.172 


-0.013 ± 0.190 


0.227 ± 0.167 


*^i<:*+{1430) 
"ii'*+(1430) 
«p0(770) 


71.373 ± 5.147 
—2 507 ± 144 
n nriQ 4- n mn 

— U.UUo ± U.UoU 

0.397 ± 0.171 
1.163 ± 0.133 


38.553 ± 4.845 
—2 983 ± 153 

U.Uol =t U.UoD 

0.636 ± 0.156 
1.394 ± 0.183 


42.430 ± 5.800 
3 294 ± 144 
n 1 1P. 4- n nof; 

0.483 ± 0.193 
1.737 ±0.145 


48.990 ± 5.035 
—2 829 ± 149 
u.zoy ± u.uyo 
0.874 ± 0.157 
1.200 ± 0.169 


V(770) 


-2.528 ± 0.306 


-3.934 ± 0.394 


1.418 ± 0.229 


-3.479 ± 0.523 


CpO(77o) 


-0.017 ±0.115 


0.083 ± 0.119 


0.035 ± 0.079 


0.057 ± 0.129 


"pO(770) 


0.350 ± 0.148 


0.398 ± 0.131 


0.406 ± 0.136 


0.457 ±0.141 


«/o(980) 


31.840 ± 2.697 


32.957 ± 2.837 


46.640 ± 3.527 


32.410 ± 2.707 


^/o(980) 


-1.309 ± 0.258 


-1.973 ± 0.284 


-2.783 ± 0.213 


-1.850 ± 0.292 


C/o(980) 


0.032 ± 0.083 


-0.000 ± 0.084 


0.006 ± 0.055 


-0.005 ± 0.085 


^/o{980) 


0.221 ±0.115 


0.258 ±0.121 


0.259 ± 0.100 


0.248 ±0.117 


a/2 (1270) 


0.101 ±0.020 


0.121 ±0.021 


0.091 ±0.022 


0.116 ±0.021 




-0.303 ± 0.235 


-0.880 ±0.217 


-0.710 ±0.311 


-0.829 ± 0.240 


CRest 


-0.042 ± 0.043 


-0.016 ± 0.043 


-0.022 ±0.037 


-0.029 ± 0.045 


finest 


n AQO _|_ n 1 74 


n fin^"^ + n 1 91 


n Qsn _|_ n 1 1 c 
u.oou in u.-Lxu 


n KCQ _|_ n 1 1 c: 


a/x(i300) 




1 '\ 748 + 1 ^Q^ 
±o. I 4:0 in x.uc/tj 




i-O.uou in i.ot I 


^/x{1300) 


n QQO + n 288 




n KVQ _|_ n 978 


n 047 _|_ n 959 


^(ii'gvr-)NR 
a(7r+7r~)NR 


109.661 ±9.777 
-3.591 ±0.129 
24.661 ±20.171 
-1.984 ±0.430 
42.197 ± 6.620 


104.821 ± 10.038 
-3.165 ±0.149 
35.387 ±9.859 
-1.531 ±0.331 
34.769 ± 9.004 


114.600 ±9.228 
3.163 ±0.155 

57.170 ± 12.660 
5.251 ±0.209 

53.560 ± 8.202 


103.400 ± 9.538 
-3.396 ±0.151 
42.490 ± 11.270 
-1.779 ±0.245 
40.960 ± 7.972 


^(7r+7r-)NR 


-0.230 ±0.184 


-0.515 ±0.197 


1.118 ±0.216 


-0.556 ± 0.242 


a 


0.147 ±0.032 


0.134 ±0.020 


0.089 ±0.025 


0.127 ±0.020 


-21og£ 


18472.5 


18465.0 


18458.5 


18465.9 



Consequently, Solutions 1 and 2 are treated as the solutions of this analysis. The 
high Kq'^ {lA'i{))TT~ fraction of Solution 1 is in agreement with some phenomenological esti- 
mates [48i] and may also be qualitatively favored by the total K — it S-wave phase shift as 
a function of m^Kir) when compared with that measured by the LASS collaboration jiol ]. 
As the likelihood difference is not found to be significant using an ensemble of pseudo- 
experiments corresponding to the luminosity of the data sample, we retain both solutions. 

To assess how well the fit to the Dalitz plot represents data as shown in Fig. [21 the time 
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TABLE IV: Summary of relative fractions where only the statistical errors are given. 



O 

Q 



K*+(892)7r- 
i<:o+(1430)^- 
p0(770)K0 
/o(980)i^o 
/2(1270)/^o 
/x (1300)^-0 

(K§7r+)NR^- 

(ii:g7r-)NRvr+ 
fvr+^-)NRi^g 



Total 



Sol. 1 Fraction (% 


Sol. 2 Fraction (% 


Sol. 3 Fraction (% 


Sol. 4 Fraction (% 


9.3 ±0.8 


9.0 ± 1.3 


8.7 ± 1.2 


9.1 ± 1.2 


61.7 ± 10.4 


17.4 ± 5.0 


20.5 ± 6.4 


28.5 ±7.4 


6.1 ± 1.5 


8.5 ±2.6 


12.8 ±2.8 


6.4 ±2.0 


14.3 ±2.7 


14.9 ±3.3 


28.9 ±5.9 


14.6 ±3.1 


2.6 ±0.9 


3.2 ± 1.2 


1.8 ±0.9 


3.1 ± 1.1 


2.3 ±0.8 


6.0 ± 1.6 


4.9 ± 1.7 


5.5 ± 1.5 


57.2 ± 11.4 


55.9 ± 13.3 


96.7 ±20.5 


58.2 ± 13.2 


2.9 ±4.7 


6.4 ±3.7 


24.1 ± 11.2 


9.8 ±5.4 


10.7 ±3.5 


7.7 ±4.1 


24.9 ± 8.4 


11.3 ±4.7 


167.1 ±0.2 


129.0 ±0.2 


223.2 ±0.3 


146.5 ±0.17 



and flavor-integrated Dalitz plot is divided into variable size bins so that each bin contains at 
least 25 events. A goodness-of-fit statistic for the multinomial distribution is then calculated 
as the sum over N bins, 

X' = -2j:n.\og^, (34) 

where nj is the number of events observed in the i-th bin and pi is the number of events 
in the i-th bin as given by the fit l^l- The distribution of this statistic is bounded by a 

distribution with — 1 degrees of freedom and one with N — k — 1 degrees of freedom 
where = 19 is the number of Dalitz plot fit parameters. The x^ /N value for the best fit is 
201.9/(189- 19 - 1). 

We measure the time-dependent CP parameters of KgTT^n' , 

Solution 1: -21og£ = 18472.5 

Acp{p'^{770)K'^s) = ±0.03l[J:i± 0.11 ±0.10, 

0f (p°(770)A:^) = (±20.0+1^ ±3.2 ±3.5)°, 

^cpip\770)K"s) = ±0.641°:^^ ±0.09 ±0.10, 

Acp{fo{9SO)K^s) = -0.06 ±0.17 ±0.07 ±0.09, 

</.f(/o(980)ir°) = (±12.7i^|±2.8±3.3)°, 

S^p{fo{980)Kl) = -0.43i°;^2± 0.09 ±0.11, 

^cp(i^*"^(892)7r") = -0.21 ±0.11 ±0.05 ±0.05, 

A0(A:*+(892)7r-) = (-0.7+^1^ ± 11.0 ± 17.6),° (35) 

Solution 2: -2 log/: = 18465.0 

Acp{p° {770) K^s) = -0.16 ±0.24 ±0.12 ±0.10, 
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FIG. 2: Signal enhanced Dalitz plot fit for — > K'^tt'^-k . When plotting a two-particle mass 
projection, we require the invariant mass of the other two two-particle combinations to be greater 
than 1.5 GeV/c^. The solid curve shows the fit result while the shaded and dashed curves show 
the qq and total background components, respectively. 



0f (p°(770)ir^) = (+22.8 ±7.5 ±3.3 ±3.5)°, 

'5cp(p°(770)/s:^) = +0.711°:^^ ±0.08 ±0.09, 

Acp{fo{9SO)K°s) = ±0.00 ±0.17 ±0.06 ±0.09, 

0f(/o(98O)ir°) = (±14.8l6^:?±2.7±3.3)°, 

S'Jp{fo{980)Kl) = -0.49l°J2± 0.08 ±0.10, 

^cp(^*^(892)7r-) = -0.20 ±0.11 ±0.05 ±0.05, 

A0(K*+(892)7r-) = (±14.6^^^:^± 11.0 ± 17.6),° (36) 

where the first error is statistical, the second is systematic and the third is the Dalitz plot 
signal model uncertainty. The At and raw asymmetry fit projections for 5° — >■ p'^{770)Kg 
and /o(980)ii'° are shown in Fig. [3], and the statistical correlation coefficients between 

the CP parameters of both solutions are given in Tables |V] and IVII The full correlation 
matrices are given in Tables llXtlXIVI 

Likelihood scans of (pf^ for both solutions are obtained by fixing 0^^^^ and redoing the fit. 
The statistical error for S^p is also determined from a likelihood scan that fixes Acp and 
(j)f^. Similarly, likelihood scans of A0 are also produced. In addition, we also perform scans 
that include the systematic and model errors by convolving the likelihood with a Gaussian 
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FIG. 3: Time-dependent Dalitz plot fit results for — > iC^vr^vr^ in (a), the p'^(770) region and 
(b), the /o(980) region using Solution 1. The top plots show the At distribution for and B^ 
tags indicated by the solid and dashed curves, respectively. These plots contain only good tags, 
0.5 < r < 1.0 and the dotted curve represents the background contribution. The bottom plots 
show the B^B^ raw asymmetry of the top plots, {N^o — Nj^o)/{N^o + Nj^o), where A^'^o (N^o) is 
the number of B^ {B^) tags in At. 



with width set to the quadratic sum of the systematic and model uncertainties. These are 
shown in Figs. HJITI 



SYSTEMATIC UNCERTAINTIES 



The systematic errors in the vertex reconstruction include uncertainties in the IP profile, 
charged track selection based on track helix errors, helix parameter corrections. At and 
vertex goodness-of-fit selection, Az bias and SVD misalignment. The parameters for w and 
Aw, resolution function, physics parameters, background shape and signal probability are 
varied by ilex. Where histograms were used, all bins are simultaneously but independently 
varied by their respective errors. The effect of CP asymmetry in the background is estimated 
by varying the CP parameters of the entire neutral B component within the physical range 
except for v'^S' which is varied within its known uncertainties. Toy MC samples 

showed some small fitting bias for CP parameters due to low statistics in each sample. 
We take this bias as a systematic uncertainty. The effect of misreconstructed signal events 
is also investigated using signal MC. This is achieved by comparing the fit result of the 
signal MC sample with another fit on the same sample, which required that the events were 
reconstructed correctly. The efficiency histogram also includes systematic uncertainties in 
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TABLE V: Correlation between CP parameters for Solution 1. 



o 



0&2 
o 



00 
05 



oCrj 

00 
OS 

CDi— i 



00 

+ 

* 



00 

+ 

* 



Acpip'i770)Kl) 


+1.00 












<(p0(770)K0) 


-0.01 


+1.00 










^cp(/o(980)KO) 


-0.29 


-0.01 


+1.00 








<(/o(980)i^O) 


-0.06 


+0.37 


+0.03 


+1.00 






^cp(^*+(892)7r-) 


+0.07 


+0.01 


+0.00 


-0.07 


+1.00 






-0.09 


+0.15 


+0.03 


+0.19 


-0.06 


+1.00 



TABLE VL Correlation between CP parameters for Solution 2. 



Ot/3 
o 



o 

it) 



00 
05 



00 



EtI 



CM 
00 
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CM 
Ol 
00 

+ 

* 

<l 



Acpip'{770)Kl) 


+1.00 












4>f{p'{770)Kl) 


-0.12 


+1.00 










^cp(/o(980)KO) 


-0.40 


-0.05 


+1.00 








< (/o(980)i^O) 


-0.13 


+0.43 


+0.14 


+1.00 






^cp(^*+(892)7r-) 


+0.02 


+0.00 


+0.01 


-0.05 


+1.00 




A^(ii:*+(892)7r-) 


+0.16 


+0.25 


-0.09 


+0.22 


+0.05 


+ 1.00 



correction factors due to tracking, Kg selection and particle identification. The systematics 
from these data to MC efficiency ratios are calculated from independent studies at Belle. 
Tag-side interference, which only affects CP eigenstates, comes from CP violation on the 
tag side {5l|, and is estimated with B —>■ 0*1^. We generate MC pseudo-experiments and 
perform an ensemble test to obtain systematic biases from tag-side interference [52]. In the 
signal model, the Dalitz plot mass, widths and other parameters are varied by their errors 
and the Blatt-Weisskopf barrier factors are included. The Dalitz plot model is the dominant 
systematic source and is quoted separately. We introduce additional resonances, K2{14:30), 
K*{1680), w(782), p°(1450), p°(1700) and /o(1710) into the signal model to estimate possible 
effects from other resonances not included in the nominal model. These resonances are 
included separately with their amplitudes and phases as free parameters and their CP 
parameters shared with the /2(1270), /x(1300) and the non-resonant component. The 
amplitudes of these additional resonances are not found to be significant. The shape of 
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FIG. 5: Likelihood scans of (t)f for pO(770)i^g (top) and 5° ^ /o(980)Kg (bottom) for 

Solution 1 (left) and Solution 2 (right). The solid (dashed) curve contains the total (statistical) 
error and the dotted box indicates the parameter range corresponding to itlo". 



the non-resonant component is empirically chosen, so different parametrizations are tested. 
They include modeling the non-resonant part with the tail of a Breit-Wigner, i?NR(s; a) = 
ia/{s + ia) and a power law whose exponent is a fit parameter, i?NR('5; a) = The fit 
differences from these alternate Dalitz plot parametrizations were summed in quadrature. 
The systematic errors for both solutions are summarized in Table IVIII and Table IVIIII As 
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FIG. 6: Likelihood scans of for B° p^{770)K'^ (top) and 5° ^ /o(980)Kg (bottom) for 
Solution 1 (left) and Solution 2 (right). The solid (dashed) curve contains the total (statistical) 
error and the dotted box indicates the parameter range corresponding to ztlo". The statistical error 
of Sqp is determined from these scans. 



the systematic uncertainty from the Dalitz plot model is the largest, it is quoted separately. 



CONCLUSION 

In summary, we perform a time-dependent Dalitz plot measurement of CP parameters 
in — s> KgTi^Ti^ decays and find two solutions that describe the data well. The first of 
these solutions may be preferred by external information from other measurements; however, 
we retain both solutions. This is Belle's first measurement of CP violation parameters in 
the {770) Kg channel and the first measurement of the CP parameters in 5° — > 

fo{9SO)Kg decays using a time-dependent Dalitz plot technique. There is currently no evi- 
dence for direct CP violation in S° ^ p°(770)K0, 5° -* /o(980)ir^ and B^ K*+{892)7c-, 
while mixing-induced CP violation in B^ p^{770)Kg and B^ fo{980)Kg decays devi- 
ates from zero by roughly 2a and is consistent with measurements in b ^ ccs transitions. 
We also measured the phase difference between B° K*^ {892)71^ and B^ K*^ {892)71^ , 
which may be used to extract 03 . 



20 




FIG. 7: Likelihood scan of Ac/) for Solution 1 (left) and Solution 2 (right). The solid (dashed) curve 
contains the total (statistical) error and the dotted box indicates the parameter range corresponding 
to ±la. 



TABLE VIL Summary of systematic uncertainties for Solution 1. 



O 

bJO 

O 



o 



o 



a. 



ftl 

QjT-H 



00 

o 



03 
00 

* 

< 



Vertex Reconstruction 

Flavor Tagging 
At Resolution Function 
Physics Parameters 
Background Model 
Signal Probability 

Fit Bias 
Misreconstruction 
Efficiency 
Tag-side Interference 
Signal Model 



0.068 
0.005 
0.022 
0.000 
0.052 
0.019 
0.013 
0.008 
0.025 
0.039 
0.045 



2.16 
0.11 
0.87 
0.04 
1.52 
0.63 
1.05 
0.38 
0.33 
0.06 
0.75 



0.045 
0.005 
0.011 
0.001 
0.020 
0.008 
0.004 
0.008 
0.008 
0.043 
0.015 



1.64 
0.20 
0.99 
0.03 
1.50 
0.54 
0.07 
0.34 
0.59 
0.03 
1.02 



0.039 
0.001 
0.013 
0.001 
0.011 
0.007 
0.002 
0.009 
0.016 
N/A 
0.016 



8.51 
0.62 
3.19 
0.12 
3.36 
3.92 
0.66 
0.31 
1.54 
N/A 
2.86 



Total 



0.112 3.17 0.069 2.79 0.050 10.99 



Amplitude Model 



0.097 3.53 0.091 3.33 0.053 17.61 
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TABLE VIII: Summary of systematic uncertainties for Solution 2. 
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O 
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O 








-s- 




1 
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Vertex Reconstruction 


0.055 


2.47 


0.031 


1.95 


0.041 


8.87 


Flavor Tagging 


0.006 


0.15 


0.005 


0.17 


0.001 


0.49 


At Resolution Function 


0.035 


0.74 


0.009 


1.00 


0.005 


3.16 


Physics Parameters 


0.000 


0.04 


0.002 


0.03 


0.000 


0.05 


Background Model 


0.078 


1.60 


0.020 


1.17 


0.014 


4.93 


Signal Probability 


0.041 


0.68 


0.012 


0.42 


0.005 


1.71 


Fit Bias 


0.013 


1.01 


0.004 


0.07 


0.002 


0.66 


Misreconstruction 


0.010 


0.35 


0.010 


0.20 


0.007 


0.19 


Efficiency 


0.017 


0.08 


0.010 


0.44 


0.028 


0.96 


Tag-side Interference 


0.039 


0.06 


0.043 


0.03 


N/A 


N/A 
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TABLE IX: Statistical correlation matrix for Solution 1. 
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TABLE X: Statistical correlation matrix for Solution 1. 
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TABLE XI: Statistical correlation matrix for Solution 1. 
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TABLE XII: Statistical correlation matrix for Solution 2. 
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TABLE XIII: Statistical correlation matrix for Solution 2. 
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+0.33 


-0.06 


-0.01 


-0.03 



TABLE XIV: Statistical correlation matrix for Solution 2. 
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1 

t= 
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1 
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O/x(1300) 


+1.00 


















^'/x(1300) 


+0.27 


+1.00 
















0(i<:g7r+)NR 


+0.42 


+0.17 


+1.00 














^(/f07r+)NR 


+0.09 


+0.34 


+0.13 


+1.00 














+0.03 


-0.08 


-0.01 


-0.24 


+1.00 












-0.10 


+0.05 


+0.04 


+0.57 


-0.07 


+1.00 








0(7r+7r-)NR 


+0.09 


+0.26 


+0.30 


-0.30 


-0.04 


-0.49 


+1.00 






^(7r+7r-)NR 


-0.25 


+0.44 


+0.03 


+0.34 


+0.10 


+0.34 


+0.15 


+1.00 




a 


+0.31 


+0.13 


+0.37 


+0.22 


+0.14 


-0.36 


+0.05 


-0.37 


+1.00 
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